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1 Introduction 

Over the last few years more and more attention has been addressed to the analysis of physi- 
cal nano-systems in order to realize intriguing applications, for instance in the field of quan- 
tum computation, and to investigate fundamental aspects of quantum mechanics. One of the 
most promising physical contexts is that of laser-driven trapped ions (for a review see, for in- 
stance HJED- 

An electromagnetic (e.m.) trap is a device which generates a suitable e.m. field confining 
a charged particle in a finite region of space. In particular, a Paul trap exploits an inhomoge- 
neous and time-dependent e.m. field which forces a charged particle to move approximately as 
a harmonic oscillator. The equilibrium point coincides with the center of the trap [3]. There- 
fore, a trapped ion is describable as a compound system made of a three-dimensional quantum 
harmonic oscillator (representing the motion of the ion center of mass inside the trap) and 
a few-level system (associated with the internal atomic state, i.e. with the 'relevant electronic 
levels' of the ion). Through the action of suitable driving laser fields, it is possible to coherently 
manipulate this system with a high degree of accuracy. In particular, the possibility of inducing 
couplings between the atomic degrees of freedom and the ion center of mass motion has been 
experimentally demonstrated in a very wide variety of settings ^[^HH]- 

In many experimental situations, only two atomic states are effectively involved in the dyna- 
mics of a ion trap. In fact, such two states are the only two effectively coupled by the laser fields 
driving the ion. Such effective couplings, for technical reasons, are usually realized through 
a third atomic level. Precisely, the two effective atomic levels - say |1) and |2) - are non- 
resonantly coupled to a third level - the auxiliary level |3) - and the respective two 'amounts of 
off-resonance' (i.e. the two detunings, that is the differences between the atomic Bohr frequencies 
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and the corresponding laser frequencies) are chosen to be equal [I]. The idea behind such a cou- 
pling scheme - the 'Raman scheme' - is that direct transitions from (and to) level |3) to (and 
from) the other two levels are forbidden by the energy conservation, while two-photon processes 
bringing from (and to) level |1) to (and from) level |2), through the auxiliary level |3), are possi- 
ble, and this realizes an effective coupling 1 1} «^ |2) . We will show that such a reasoning is correct 
within a certain approximation which will be clarified in the following (see Sections H] and 0) . 

In the present paper, we analyze the Raman coupling schemes realized by laser-driven trapped 
ions deducing, by means of a rigorous perturbative approach, the expressions of the relevant 
effective couplings. To this aim, we exploit a recently introduced perturbative method based on a 
suitable decomposition of the time evolution operator associated with a quantum Hamiltonian [5] 
(for the case of a time-independent Hamiltonian, see also [0113 IE])- In particular, we investigate 
the problem of determining the effective Hamiltonian - i.e. the Hamiltonian describing the 
effective couplings - in the case where two or more Raman coupling schemes are simultaneously 
active, so providing a rigorous proof of the additivity of the effective couplings. 

We will show, moreover, that in the special case of a single Raman coupling our result 
coincides with the result previously obtained by means of a time-independent perturbative app- 
roach 

The paper is organized as follows. In Section |2 we describe the general form of the Hamilto- 
nian associated with a trapped ion Raman scheme. In Section |31 we introduce the perturbative 
method which is the main tool of the paper. We apply this method in Section^ where we ana- 
lyze the dynamics of a double A Raman scheme. Finally, in Section some conclusive remarks 
are drawn. 

2 The physical system 

The general form of the quantum Hamiltonian of a trapped three-level ion addressed by a set 
of laser beams coupling the atomic level |3) with the other two levels (with suitable ion-laser 
detunings) is the following: 

H(t) = H + H B + H R (t), (1) 

where 



Hp = hwiau, H B = hv ala 



Z=l,2,3 a=x,y,z 



H R (t) = (hn 13 (t)a 13 + h.c.) + {hh 23 (t)a 23 + h.c), (2) 

with d[ m = \l) (m\, l,m = 1,2,3, {|0}f=i being the considered three atomic levels and {h\vi}f =1 
the corresponding energies ^ uj r , for / ^ r); {a a '■ a = x,y, z} are the annihilation operators 
associated with the center of mass harmonic motion along the axes x, y, z, namely, 

^ = © 1/2 (" + i^)' a * = (i) V2 ( l+ ^)' (3) 

which, of course, satisfy the well known bosonic commutation relations [do,^] = [d^djj] = 0, 

[ao-,^] = 5 aj 3, for all a,f3 = x,y,z. Without loss of generality, the three harmonic oscillator 
frequencies have been taken to be equal: v x = v y = v z = v > 0; thus, we deal with a spherically 
symmetric trap. 

We have denoted by i i — > Qj 3 (t), j = 1,2, operator-valued functions acting in the bosonic 
Hilbert space (i.e. the Hilbert space associated with the vibrational degrees of freedom); their 
specific structure is determined by the specific laser configuration. For instance, for a single 
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'A Raman coupling' - involving two laser beams with complex strengths (proportional to the 
laser amplitudes and to the atomic dipole operators, and including the laser phases), wave vectors 
and frequencies 513, ki 3 , ui 3 and g 23 , k 23 , ^23> respectively - with the lasers in a 'travelling wave 
configuration', we have: 

Ai3(t) = n 9l3 e- i (^^ t ), Cl 23 (t) = hg 23 e<^- r ^ t ), (4) 

where the laser frequencies are fixed in such a way that the two couplings share the same 
ion-laser detuning 

A = UJ 3 - LOl - U) 13 = UJ 3 - UJ 2 - UJ 23 ^ °) 

in order to allow the typical 'two-photon processes' |1) <^ |2) of the Raman scheme. In (jlj the 
vector operator r := (x,y,z) is the ion center of mass position operator and its presence in the 
interaction Hamiltonian H^(t) is responsible for the interaction between atomic and vibrational 
degrees of freedom. The link between this operator and the annihilation and creation operators 
stems from relations (j2J). With regard to the coefficients 513 and g 23 appearing in the definition of 
the operators &i 3 (t) and ^23^) respectively, we recall that they are given by gj 3 := —h~ 1 dj 3 -Ej 3 , 
where Ej 3 is the complex amplitude (i.e. including information about the initial phase of the 
field) of the laser field tuned near the \j) — > |3) Bohr frequency, while dj 3 is the atomic dipole 
operator matrix element involving the atomic states \ j) and |3). 

A relevant feature of the laser configuration specified by relations (@J) is that, in this particular 
case, the Hamiltonian H (t) can be transformed, by passing to a suitable rotating frame (i.e. by 
switching to an ad hoc interaction picture), into a time-independent Hamiltonian, which can be 
then treated by means of a time-independent perturbative approach, see 

However, more complicated laser configurations are possible and useful for various applica- 
tions, and, in general, one cannot find a simple rotating frame where the Hamiltonian of the 
system becomes time-independent. Therefore, it will be convenient to apply a time- dependent 
perturbative approach. 

3 The time-dependent perturbative approach 

In order to study the class of dynamical problems associated with a quantum Hamiltonian of the 
form (^Q), one can fruitfully exploit a time-dependent perturbative method based on a suitable 
decomposition of the evolution operator which is a generalization of the classical Magnus 
expansion [ID], 

Consider a quantum system whose Hamiltonian is made of two components, the unperturbed 
energy operator Hq, and a perturbation H (X;t), in general time-dependent: 

H(X;t) = H + H <> (X;t). 

We assume that A 1— > H (X;t) is an analytic function of the (real) perturbative parameter X 
(#o(0;t) =0). 

Introducing the evolution operator generated by the unperturbed component Hq - namely, 
Uq := e~n Hot - and the the evolution operator T(X;t) associated with the interaction picture 
Hamiltonian H(X;t) := Uo(ty H<>(X; t) Uo(t), it is possible to factorize the total evolution opera- 
tor U(X; t) of the system as 



U(X;t) = U (t)T(X;t). 



(5) 
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We can now consider the following exact decomposition of the interaction picture evolution 
operator T(A; t) as a product of unitary operators: 

T(A; *) = exp (-i Z(X; t)) exp f-i C(A; t) dt^ exp (i Z{X)) , (6) 

where Z(X;t), C(X;t), Z(X) are selfadjoint operators and 

Z(X) = Z(X; 0), Z(0; t) = C(0; t) = 0, Vt 

It is worth noting that no time-ordering operator is present in (JBJ. Moreover, observe that in 
the special case where Z(X;t) = 0, Vt, we have the Magnus expansion, provided that C(A;t) is 
regarded as the time derivative of the Magnus unitary operator generator, see ^U]; but, as we 
will see below, decomposition (jSJ is actually a generalization of the Magnus expansion. 



3.1 Imposing a gauge condition 

It can be shown that the operators C(X;t) and Z(X;t) are not uniquely determined in decom- 
position @: in fact, there are infinite possible solutions, namely, solutions compatible with the 
general form of such decomposition. A unique solution can be singled out by imposing an addi- 
tional 'gauge condition'. A typical example is the case where the interaction picture Hamiltonian 
H(X; t) is a almost periodic 1 operator-valued function of time, in particular, an operator-valued 
trigonometric polynomial with respect to the time variable. In this case, a remarkable gauge is 
fixed by the following tern of conditions (see 

1) C(X;t) = C(A;0) = C(A), V t; 

2) the function t i— > Z{X\ t) satisfies: 

lim t^Zfct) = 0; 

t— >oo 



3) the mean value of the function 1 1— > Z(X;t) is zero, namely: 

lim - I Z(X;t)dt = 0. 
r^oo r J 

As it will be seen in Section^ this is precisely the case occurring in our applications, due to 
the fact that the operators ^13 (t), &23(t) are indeed trigonometric polynomials: 

Kj 

£l j3 (t) = J2fy 3 e i "" t , Rj £ N, {cjj, . . . , ujj 3 } CR, j = l,2, 

K=l 

where {fi},...,fiW = are operators acting on the bosonic degrees of freedom; see, for instance, 
relations @. 

3.2 Determination of C(A;t) and Z(X; t) up to a gauge condition 

Using the formula - reported, for instance, in ^2] ~ f° r the derivative of the exponential of an 
operator- valued function 1 1— ► F(t), i.e. 

X A standard reference on almost periodic functions is 
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we can write the Schrodinger equation for the interaction picture evolution operator as 

H{X;t)T{X;t) = ihf(X;t) = he~ iZ ^ (e isZ{x > t] Z{\- t)e~ isZ{x ^ ds e -^o c(A;t)d V z(A) 

+ he -iZ^t)j^sf<C(X:i) dt c(A . t)e w/*0(A;i) dt) dg ^if^i) dt^A) (?) 

Applying to each member of equation Q the operator e J ^( A; *) on the left and the operator 
e -iZ(\) e if c(\;t) dt Qn ^ e right, we get the following equation relating the operators C(A;i) and 
Z{X\ t) with the interaction picture Hamiltonian: 

Ad cxp(iZ(A;t)) ^(A;t) = h (Ad cxp(isZ(A;t)) Z(A;f) + Ad cxp ^ is ^ c{A;t)d ^C(A; t)j ds. (8) 

We recall that, given linear operators X (invertible) and Y , Ad^Y := XYX~ . If the operator X 
is of the form X = e x , we can use the well known relation 

oo ^ 

Ad cxp(x) Y = exp(ad x )Y = £ ^d^Y, (9) 

m=0 

with ad™ denoting the m-th power (ad^- = Id) of the adjoint super-operator adx defined by 
adx Y := [X, Y]. Applying formula © to equation ©, and performing the integrals, we obtain: 

OO '7TJ, '7TI 

V l — ad? . H(X; t) = h y y ad?, + , Z(X; t) 

m \ Z(X;t) ^ ' > ( m+ i)| Z(A;t) V ' ' 

m=0 m=0 v ' 

OO / -\ 772 

m=0 

Next, using the Taylor expansions of the operator- valued functions A i— ► H (A; i), A h- > Z(A; i) 
and A i-» C(A;t) (recall that H(0;t) = Z(0;t) = C(0;t) = 0), i.e. 

OO OO OO 

H(X;t) = J2* n Hn(t), Z(X;t) = J2^ z n(t), C(A; i) = £ X n C n (t), (11) 

n=l n=l n=l 

from formula (|l(Jj) we can determine (in a non-unique way) - order by order with respect to the 
perturbative parameter A - the operators {Z n (t)}^° and {C n (i)}f°. 

3.3 Solution corresponding to conditions [TJHHI) 

As anticipated, decomposition © - or, equivalently, equation (|1U|) - admits infinite solutions. 
Therefore, it is necessary to impose a gauge condition in order to obtain a specific solution. For 
instance, as already mentioned, the condition Z(X; t) = 0, Vi, allows to obtain a precise solution, 
that is the well known Magnus expansion. 

In the following, we will use, instead, the gauge fixed by the tern of conditions ^)^HJ)- With 
these conditions, one can single out a unique solution for the operators C(X;t) and Z(X;t); 
precisely, we get an infinite set of equations that can be solved recursively, order by order, for 
obtaining the operators {C n (t) = C n (0) =: C n }^ =1 and {Z n (t)}^ =1 that appear in the power 
expansions (|lljl. up to an arbitrary order N; specifically, it turns out that 

1 st order: 

Ci = lim - f JT^HiMdt, 

T ^°° T Jo 
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Zx(t) = J h- l H x {i)&i - dt- Urn ]- J (j (h^H^t) -Ci) dt\ dt; (12) 
2 nd order: 

C 2 = Hm i jT r- ad Zl(t) (V^t) + d) + h^H^A dt, 
Z 2 (t) = jf f 1 ad Zl(i) (h^Hiit) + Ci) + ft" 1 ^)) dt - C 2 t 

- Hm i jT Qf* Q ad Zl(t) (V^t) + d) + fc-^t) - C 2 ) dt) dt; (13) 



Notice that the operators C\, Z±(t), C 2 , Z 2 (t), . . . are indeed selfadjoint, coherently with our 
previous assumption. 

As already mentioned, a typical case occurring in applications (and, in particular, in the 
application considered in the present paper) is the case where the coefficients of the perturbative 
expansion of the interaction picture Hamiltonian are operator-valued trigonometric polynomials 
with respect to the time variable (hence, almost periodic functions of time). In this case, the 
gauge fixed by conditions is such that the functions {t h- > Z n (t)}^ =1 are zero- mean- valued 
trigonometric polynomials; as a consequence, all the 'secular terms' are concentrated in the 
component of the perturbative decomposition of the evolution operator which is generated by 
the operators {C n }^ =1 , i.e. in the one-parameter group of unitary operators {exp(— iC(X)t)}t£R- 



3.4 N-th order truncation of the perturbative decomposition (Jf)]) 

Once that the operators C\, Zi(t), . . . have been obtained recursively up to a certain perturbative 
order N > 1, one can write the following iV-th order approximation of the interaction picture 
evolution operator: 

T(X;t) « exp (-iZ (N) {\; t)) exp (-iC (N) (A)t) exp (iZ m {X)) , (14) 

with 

N N 

C W (A) := A n C n and Z (JV) (A;t) := X n Z n (t). 

n=l n=l 

We stress that the iV-th order truncation (|14(1 preserves the fundamental unitary nature of the 
interaction picture evolution operator T(A;t). 

From formula (|14|) we find that the overall evolution operator of the system admits the 
following iV-th order approximation (recall relation ©): 

U(X; t) « U (t) exp (-iZ (JV) (A; t)) exp(-iC (JV) (A) t) exp (iZ (N) {X)) 

= exp ( - iZ (JV) (A; t)) U (t) exp {-iC iN) (X)t) exp (iZ (JV) (A)) , (15) 

where Z (JV) (A;t) = U {t)Z m (X; £)C/ (*) t - From relation O it follows that 
£/(A;t) ^exp(-iZ {N \X;t))U^\X;t) exp (iZ (JV) (A)) , 

where U^\X;t) = Uo(t) exp (— iC (JV) (A)t) can be regarded as the evolution operator associated 
with the effective Hamiltonian 

H ( e £(X;t) = H + h& N) (X;t) with C m (X;t) := oi C (Ar) (A)e^>*. 
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Thus, at the iV-th perturbative order, N > 1, the total evolution of the system can be 
decomposed into a 'dynamical dressing' - i.e. the passage to a generalized interaction picture 
generated by the time-dependent transformation exp ( - iZ (JV) (A; t)) , see - and the evolut ion 
generated by an effective Hamiltonian H^P (A;i) having the fundamental property that the 
corresponding interaction picture Hamiltonian, with respect to the reference Hamiltonian Hq - 
namely, C (JV) (A) - is time-independent. 

4 Raman schemes: effective coupling 

The perturbative approach described in the preceding section turns out to be a powerful tool for 
studying Raman schemes and for deducing the 'effective couplings' |1) <^ |2). For the sake of 
definiteness, we consider the situation where two Raman setups are simultaneously present. This 
case contains that of a single Raman coupling as a special case 2 . Therefore, a direct comparison 
of the behaviors associated with one or two Raman schemes can be given. It is worth noting 
that the results obtained for two Raman schemes may be immediately generalized to the case 
of several Raman couplings. Therefore, the situation under scrutiny allows to illustrate all the 
relevant conceptually remarkable aspects without introducing any cumbersome notation. 

The Schrodinger picture Hamiltonian describing two simultaneously active Raman schemes 
is given by equations Q and ©, now taking 

%3 W = g j3 e-<^-r-^t) + g y i(k' j3 -r- J : n j = ^ % 

where gj 3 , ujj 3 , kj 3 , with j = 1,2, are the coupling constants, frequencies and wave vectors 
associated with the first couple of Raman lasers, and g'j 3 , u^ 3 , kj 3 the analogous quantities for 

the second Raman scheme. We can assume that gj 3 7^ (and, of course, ujj 3 7^ 0, kj 3 7^ 0), 
j = 1, 2, so that the special case of a single Raman coupling is recovered for g'j 3 = 0. As already 
mentioned, the operator f is the ion center of mass position operator. 

In order to generate 'two-photon processes' involving levels |1) and |2) (which is the main 
feature of the Raman coupling), the ion-laser detunings are fixed in such a way that 

(u 3 - UJi) - CJ13 = (uj 3 - L0 2 ) - W23 = A / 0, 
(cj 3 - ui) - uj' 13 = (uj 3 - uj 2 ) - ^23 = A' 7^ 0. 

In the following, we will consider the regime where A 7^ A'; hence, the special case of a single 
Raman coupling is recovered only for g'-^ = 0, j = 1,2. We can also assume, without loss of 

generality, that |A| > |A'|, for g'- 3 ^ (and, of course, ujj 3 7^ 0, 7^ 0), j = 1,2. 
We will further suppose that the following high detuning conditions are satisfied: 

|A| » hsUtoU (|A| >) |A'| » biaUskl. ( 16 ) 

Passing to the interaction picture with respect to the reference Hamiltonian Hq, the Schro- 
dinger picture Hamiltonian H(t) is transformed into the interaction picture Hamiltonian 

H(t)=H B + H R (t), 

where 

H R {t) = (nn 13 (t)a 13 + h.c.) + (hn 23 {t)& 23 + h.c), 

2 In Section El we have mentioned that the case of a single Raman scheme can be exceptionally treated by 
a time-independent perturbative method. 
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1,2. 



At this point, in order to put in evidence the 'natural perturbative parameter' of this model, it 
will be convenient to introduce the dimensionless (interaction picture) Hamiltonian Sj(X;t) by 
setting: 

(hAr'Hit) =:Sj(X;t) 

= A £ xaia a + A £ (x^"^^^ + x^-^V^)^ + h.c.) , 

a=x,y,z j=l,2 

where A is the dimensionless (real) perturbative parameter defined by 

A := g = max{i/, |#i 3 |, [523 1, |5i 3 |, I023IK 

and k = x i3 = g j3 /g, k[ j3 = g' j3 /g, j = 1, 2. 

We notice explicitly that < x < 1, |xj 3 | < 1, |x^ 3 | < 1, j = 1,2, and, due to conditions ((To]) , 
we have that |A| <C 1; hence, A is indeed a 'good perturbative parameter' (compare with the 
time-independent perturbative approach used in |^ for the single Raman coupling). 

Next, applying formulae (j!2j) and l)13J) to the interaction picture Hamiltonian /iAf}(A;t), we 
easily find the following expressions for the operators C\, C 2 (which, as we have seen, determine 
the second order effective Hamiltonian): 

XC 1 = h~ 1 H B = v (17) 

a=x,y,z 



A^C 2 = CJiffn + W 2 <T 2 2 + £> 3 <7 33 + ( ( 9ue %KW ' T + 01 2 e ifc i2 r ) <T12 + h.c. ) , (18) 

where we have set 



^■•="^ J "1.2, W 3 .= £ + X , : 

and 



513532 / 5i 3 532 T T T Ti T> fi rm\ 

512 := A , 5i2 := A , » fc i2 := % - ^23, k 12 := k 13 - k 23 , (19) 

with g 3j = g* 3 , g' 3j = g'* 3 , j = 1,2. 

Now, according to what we have shown in Subsection 13.41 the complete temporal evolution of 
the Raman-driven trapped ion, evaluated within the second order in the perturbative parameter, 
can be written as 

U(X,t) « e-^ 0i exp(-iZ (2) (A;t))e- i ( ACl+A2c2 ) t exp(iZ (2) (A)) 

= exp(-ii (2) (A;t))e-t^* e -*(AC7 1 +A2c 2 )t exp( - Z (2) (A))i (2Q) 

where Z (2) (A; t) := e~n Hot Z^(X; t)eH Hot is an operator- valued function of time, which, as a con- 
sequence of of conditions , is a zero- mean- valued trigonometric polynomial, hence contains 
only 'oscillatory terms'; its explicit form, however, is not relevant for the purposes of the present 
paper and will be omitted. 

Thus, as relation (|2T)|) shows, up to the 'dynamical dressing' generated by the time-dependent 
unitary transformation exp(— iZ(X; t)) ~ exp(— iZ (2) (A;t)), the temporal evolution of the system 
is described, at the second perturbative order, by the evolution operator e -ft^ot e -«(ACi+A 2 C2)^ 
Such an evolution operator may be regarded as a Schrodinger picture evolutor factorized into the 
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the unperturbed evolution e~n Hot and the interaction picture evolution e - i (^Ci+X 2 C 2 )t ^ 

is the one-parameter group of unitary operators generated by the time-independent 'interaction 

picture Hamiltonian' \ C\ + X 2 C 2 . 

The corresponding Schrodinger picture Hamiltonian is then given by 

:= H + e-^o* (Ad + X 2 C 2 ) ei B « 

and - as it is easily checked using formulae (|17|) and ()18|l - it can be decomposed by means of the 
complementary orthogonal projectors Pi 2 = <5"n + cr 22 and 033 (-P12 + ^33 = Id, [-Pi2 5 o"33] = 0): 

H%(t)=H$(t)^ + H%(t)d Sa , [H$(t),P 12 ] = [H^(t),a 33 ] =0, 

where 

^df (*)A2 = ^ X] "1^2 + M w l + ^l)* 5 "!! + ^(^2 + ^2)<5-22 

+ ((^lae^^ 12 -^ 12 *) + %i 2 e-<^2-^i 2 t)^ ^ + h c ) j (21) 

H { e ft{t)a 33 = hv ^2 ala a °33 + ^3 + w 3 )<T3 3 , 
a=a;,2/,z 

with the effective frequencies u±2, u/ 12 defined by 

U12 := U 2 -UJl = L> 13 - U 2 3, v'12 ■= U)' 2 ~ U>[ = Uj[ 3 - LV 23 - (22) 

The result found has a simple and transparent interpretation. In fact, the effective Hamil- 
tonian H ( ^(t) is the sum of two completely decoupled Hamiltonians, 'living' respectively in the 
ranges of the orthogonal projectors P\ 2 and a 33 . It is worth noting the remarkable fact that the 
Hamiltonian (t) P\ 2 can be regarded as the standard Hamiltonian of a trapped two-level ion 
in interaction with a couple of laser fields characterized, respectively, by the following effective 
parameters: g\ 2 , k± 2 , oj\ 2 and g' 12 , k' 12 , u[ 2 (see formulae (|T9]l and (|22|)). 

5 Conclusions 

We can eventually draw the following conclusions: 

• the quantum Hamiltonian describing the physical system of a trapped three-level ion in 
interaction with a set of laser beams generating Raman couplings can be successfully 
treated by means of a suitable time-dependent perturbative approach; 

• the result of this treatment, in the special case of a single Raman coupling, coincides with 
the result obtainable by means of a time-independent perturbative approach, see 

• at the second perturbative order (hence, with a high degree of accuracy), the dynamics 
of the system is given, up to a 'dynamical dressing', by the dynamics associated with 
an effective Hamiltonian describing two decoupled subsystems: an effective laser-driven 
trapped two-level ion and a simple harmonic oscillator; this result shows also the additivity 
of the effective couplings: indeed, the effective Hamiltonian associated with a certain set 
of Raman couplings is characterized by a coupling term which is the sum of the single 
effective couplings (see formula (|21j)): 

• although a detailed analysis of this point is beyond the aim of the present paper, it is 
worth mentioning the fact that the behavior of certain experimentally observable quanti- 
ties, in correspondence to certain initial conditions of the system, is rather well reproduced 
by the effective dynamics only - i.e. neglecting the effect of the dynamical dressing - as 
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a consequence of the temporal coarse-graining introduced by the experimental apparatus 
(which has, unavoidably, a fixed temporal resolution); nevertheless, due to the presence of 
the dynamical dressing, in addition to the standard Rabi oscillations between the effective 
levels |1) and |2), also fast transitions coupling these two levels with the auxiliary level |3) 
take place; then, if the auxiliary level is an excited level with non-negligible decay rate to- 
wards (at least one of the) levels |1) and |2), the fast transitions to the auxiliary level, com- 
posed with decays, gradually injects decoherence into the effective coherent cycle involving 
levels |1) and |2) (see [2]), result which is actually in agreement both with experimental 
observations |131 114j and with numerical simulations on the basis of a phenomenological 
master equation [IS] . 

A detailed study of the experimentally observable implications of the dynamical dressing, 
especially with regard to the decoherence effects, is still work in progress ^j, and it seems to 
raise very intriguing issues both on the theoretical and on the experimental side. 
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